Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



An unusual Ha nebula around the nearby neutron star 

RX J1856.5-3754* 

M. H. van Kerkwijk^ and S. R. Kulkarni^ 

^ Astronomical Institute, Utrecht University, P. O. Box 80000, 3508 TA Utrecht, The Netherlands 
^ Palomar Observatory, California Institute of Technology 105-24, Pasadena, CA 91125, USA 

Submitted 2 September 2001 

Abstract. We present spectroscopy and Ha imaging of a faint nebula surrounding the X-ray bright, nearby 
neutron star RX J1856.5— 3754. The nebula shows no strong lines other than the Hydrogen Balmer lines and 
has a cometary-like morphology, with the apex being approximately 1" ahead of the neutron star, and the tail 
extending up to at least 25" behind it. We find that the current observations can be satisfactorily accounted for 
by two different models. In the first, the nebula is similar to "Balmer-dominated" cometary nebulae seen around 
several radio pulsars, and is due to a bow shock in the ambient gas arising from the supersonic motion of a 
neutron star with a relativistic wind. In this case, the emission arises from shocked ambient gas; we find that the 
observations require an ambient neutral Hydrogen number density rijjo — 0.8 cm~^ and a rotational energy loss 
_E ~ 6 X 10"^^ ergs~^. In the second model, the nebula is an ionisation nebula, but of a type not observed before 
(though expected to exist), in which the ionisation and heating are very rapid compared to recombination and 
cooling. Because of the hard ionising photons, the plasma is heated up to ~ 70000 K and the emission is dominated 
by coUisional excitation. The cometary morphology arises arises naturally as a consequence of the lack of emission 
from the plasma near and behind the neutron star (which is ionised completely) and of thermal expansion. We 
confirm this using a detailed hydrodynamical simulation. We find that to reproduce the observations for this case, 
the neutral Hydrogen number density should be njjo ~ 3 cm~^ and the extreme ultraviolet flux of the neutron star 
should be slightly in excess, by a factor ~ 1.7, over what is expected from a black-body fit to the optical and X-ray 
fluxes of the source. For this case, the rotational energy loss is less than 2 x 10^^ ergs~^. Independent of the model, 
we find that RX J1856.5— 3754 is not kept hot by accretion. If it is young and cooling, the lack of pulsations at 
X-ray wavelengths is puzzling. Using phenomenological arguments, we suggest that RX J1856.5— 3754 may have 
a relatively weak, few 10^^ G, magnetic field. If so, it would be ironic that the two brightest nearby neutron stars, 
RX J1856.5— 3754 and RX J0720.4— 3125, may well represent the extreme ends of the neutron star magnetic field 
distribution, one a weak field neutron star and another a magnetar. 
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1. Introduction 

From studies with the EINSTEIN X-ray satellite and 
the ROSAT all-sky survey, nearly a dozen bright soft 
X-ray sources have been i dentified wit h nearby neutron 
stars ( ICaraveo et al. 1996| ; [Motch 2000| ). The hst includes 
nearby ordinary pulsars (e.g., PSR B0656-I-14), a millisec- 
ond pulsar (PSR J0437— 4715) and a gamma-ray pulsar 
(Geminga, likely a radio pulsar whose beam does not in- 
tersect our line of sight). Apart from these, the sample 
also includes six radio-quiet, soft X-ray emitting neutron 
stars whose origin is mysterious ( Treves et al. 200"o| ). 
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Observatory, Paranal, Chile (ESO Programmes 63.H-0416 and 
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The emission from all six sources appears to be com- 
pletely thermal. This might make it possible to derive ac- 
curate temperatures, surface gravities, and gravitational 
redshifts from fits to atmospheric models. Therefore, there 
has been growing interest in the using these neutron stars 
as natural laboratories for the physics of dense matter. 
In particular, Lattimcr & Prakash (2001) show that a ra- 
dius measurement good to 1 km accuracy is sufficient to 
usefully constrain the equation of state of dense matter. 



Of the six sources, the brightest is RX J1856.5— 3754 
( [Walter et al. 1996| ). [Walter fc Matthews (1997D 



used the Hubble Space Telescope {HST) to i dentify 
RX J1856.5— 3754 with a faint blue star, for which Waltei 
measured a distance of about 60 pc by direct 



(20011) 



trigonometric parallax. [Pons et al" (200l| ) showed that the 
totality of the data (optical, UV, EUV and X-ray) can be 



2 



M. H. van Kerkwijk and S. R. Kulkarni; An unusual Ha nebula around RX J1856.5— 3754 



accounted for by a neutron star with a surface tempera- 
ture, fcTcff — 50 eV and radius ~7 (d/60pc) km. 

As yet, no pu lsations have been detected for 
RX J1856.5-3754 (^ons et al. 200l|; iBurwitz et al 



200l[)~| In contras t, for the second brightest source 
RX J0720.4-3125, [Haberl et al. (1997D ) found clear pul- 
sations, with a period of 8.39 s. This source has also been 



identified with a faint blue star (Motch & Haberl 1998; 
Kulkarni & van Kerkwijk 1998). The long rotation period 
makes the nature of RX J0720.4— 3125 rather mysterious. 
Models that have received some considerations are a neu- 
tron star accreting from the interstellar medium and a 
middle-aged magnetar powered by its decaying magnetic 
field (ibid.; |Heyl fc Kulkarni 19981 ). 

In order to be able to secure accurate physical pa- 
rameters from observations of these sources, we need 
to understand their nature. In particular, we need to 
know the composition of the surface and the strength 
of the magnetic field. Motivated thus, we have em- 
barked on a comprehensive observational programme on 
RX J1856.5— 3754. In our first paper (Van Kerwijk & 



Kulkarai 2001, hereafter Paper I), we presented the opti- 
cal spectrum and optical and ultraviolet photometry from 
data obtained with the Very Large Telescope (VLT) and 
from the HST archive. We found that all measurements 
were remarkably consistent with a Rayleigh- Jeans tail, 
and placed severe limits to emission from other mech- 
anisms such as non-thermal emission from a magneto- 
sphere. 

In this paper, the second in our series, we report upon 
a cometary nebula around RX J1856.5— 3754, which was 
discovered during the photometric and spectroscopic ob- 
servations reported on in Paper I. We find that the neb- 
ula could be a pulsar bow shock, similar to what is seen 
around other pulsars, or a photo-ionisation nebula, of a 
type that has not yet been seen anywhere else, but which 
arises naturally in mode ls of neutron star s accreting from 
the interstellar medium ( Blaes et al. 1995 ). For either case, 
the nebula offers a number of new diagnostics: the density 
of the ambient gas and the three dimensional motion of 
the neutron star, and either a measurement of the rota- 
tional energy loss E from the neutron star, which so far 
has not been possible to make for lack of pulsations, or a 
measurement of the source's brightness in the extreme ul- 
traviolet, which is impossible to obtain otherwise because 
of interstellar extinction. 

The organisation of the paper is as follows. First, in 
Sect. [2], we briefly summarise the relevant observed prop- 
erties of RX J1856.5— 3754. Next, in Sect. ^, we present 
the details of the observations leading to the discovery of 
the nebula, and in Sects. ^ and || we investigate its nature. 

Given the detailed modelling that we have undertaken, 
the article is necessarily long. Sect. || offers a summary of 
our results and contains a discussion of the clues about the 
nature of this enigmatic source offered by the discovery 
of the nebula. We recommend that those of less patient 
inclination read Sect. ^ first. 



2. RX J1856.5-3754 

The spectral energy distribution of RX J1856.5— 3754 has 
been studied in detail by Pons et al. (200l[ ), using opti- 
cal and ultraviolet measurements from HST and soft X- 
ray measurements from EUVE and ROSAT. They find 
that the overall energy distribution can be reproduced 
reasonably well with black-body emission with a temper- 
ature kToo = 48 ± 2eV, tt-h = 2.2^°;^ x lO^Ocm-^ and 
Roo/d = O.lliO.Ol kmpc^"'^. The black-body model is not 
a very good fit to the X-ray data, which require a some- 
what higher temperature and lower column (63±3eV and 
1.0 ± 0.2 X lO^'^cm^, respectively; [Burwitz et al. 2001| ). 
This indicates that the real atmosphere is more compli- 



cated. Indeed, Pons et al. (2001) find that more detailed 
atmospheric models for compositions of 'Si-ash' and pure 
Fe give better fits. For those, however, one expects strong 
lines and bands, which appear not to be present in high- 
resolution X-ray spectra obtained with Chandra ( Burwitz 
let al. 2001| ). 

In Paper I, we reported VLT photometric and spec- 
troscopic observations of RX J1856. 5-3754. Our VLT 
spectrum, over the range 4000-7000 A, did not show any 
strong emission or absorption features. With considerable 
care to photometric calibration, we obtained photometry 
from VLT and archival iJST" imaging data. We found that 
over the entire optical through ultraviolet range (1500 A- 
7000 A), the spectral energy distribution was remarkably 
well described by a Rayleigh-Jeans tail, f\ cx A^'' with 
Ay = 0.12±0.05 and an unextincted spectral flux of /aq = 
(3.36±0.17) X 10-19 erg s'^ cm'^ A"^ at Aq = 5000 A. The 



reddening corresponds to Nh — 2.4 x 10^" cm^^ (Prcdchl 
|fc Schmitt 1995 ), consistent with that estimated from X- 
ray spectroscopy (see above). 

Walter (2001) used HST to measure the proper mo- 
tion and parallax of RX J1856.5— 3754; he found /x = 
332 ± Imasyr-i at position angle 9 = 100?3 ± 0?1 and 
TT — 16.5 ± 2.3 mas. The parallax corresponds to a dis- 
tance of 6llg pc. The proper motion is directed away from 
the Galactic plane, (^;, ^f,) = (63.8, —325.8) m asyr"^ and 
away from the Upper Scorpius OB association. Walter sug- 
gests that RX J1856.5— 3754 may have been the compan- 
ion of C Oph, which is a runaway star that likely origi- 
nated in Upper Sco, but Hoogerwerf ct al. (2001 ) argue 
that PSR J1932-hl059 was the companion to C Oph. In 
any case, the proper motion places RX J1856.5— 3754 in 
the vicinity of Upper Sco association about a million years 
ago. Thus, it is quite plausible that RX J1856.5— 3754 was 
born in the Upper Sco association, in which case its age 
is about a million years. 

The radial velocity of RX J1856.5— 3754 is unknown. 
However, assuming an origin in the Upper Sco association, 
it should be about -55kms"i ( [Walter 200l| ). If so, the 
angle between the plane of the sky and the velocity vector 
is « ~ 60°. However, the radial velocity (and thus i) does 
depend on the distance of RX J1856.5— 3754. For instance, 
if the distance is 100 pc, then the radial velocity required 
for an origin in Upper Sco is about — 20kms~i. 
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Fig. 1. Spectrum of the nebula around RX J1856.5— 3754, 
showing emission in Ha and H/3, but not in any other 
lines. The continuum emission from RX J1856.5— 3754 has 
been subtracted. The slight rise in the continuum to longer 
wavelengths is due to a faint, red extended object close to 
RX J1856. 5-3754 (see Sect. and Fig. |). 



3. Observations 

In 1999 June, we used the VLT to obtain optical spectra 
of RX J1856.5— 3754. The reduction of these data and the 
resulting spectra are described in detail in Paper I. In 
the spectra, we discovered extended Ha and H/3 emission 
around the source. We formed spectral images along the 
two slits which we used for the observations ~ one roughly 
along the proper motion direction and another roughly 
perpendicular to it (over stars F and L, respectively; see 
Fig. ^ and Fig. 1 in Paper I). From these images, it was 
clear that the nebula is most extended along the path that 
the object has taken. 



3.2. Ha image 

To study the nebula in more detail, follow-up Ha imag- 
ing was done for us in Service time with the VLT in May 
2000. In addition, deep B and R-band images were ob- 
tained. The reduction of these images and the B and R- 
band photometry are described in Paper I. The reduced 
Ha image is shown in Fig. ^ In order to bring out dif- 
fuse emission better, we also show an image in which the 
stars have been removed by point -spread function fitting 
(which was produced in the process of doing photometry 
with DAOPHOT; see below). Apart from the faint Ha neb- 
ula (inside the box), also more extended diffuse emission 
is present. The latter also appears in our B- and R-band 
images and hence does not reflect Ha emission, but more 
likely scattered light, perhaps due to dust associated with 
the R CrA complex. 

Figure also shows enlargements around 

RX J1856. 5-3754. While the nebula is clearly visi- 
ble in the Ha image, it is hard to judge its extent because 
of the crowded field. Furthermore, a faint star and 
small extended object distort the shape of the nebula's 
head (the latter also contributes to the long-wavelength 
continuum of the nebular spectrum shown in Fig. |^). 

We tried to remove the contribution of broad-band 
sources by subtracting a suitably scaled version of the R- 
band image. For this purpose, we used the image differ- 
encing technique of Alard fc Lupton (199^ ), in which the 
better-seeing image of a pair (in our case, the R-band im- 
age) is optimally matched to the worse-seeing image (Ha) 
before the difference is taken. The matching is done by 
constructing a kernel, convolved with which the PSF of 
the better-seeing image matches the PSF of worse-seeing 
image optimally (in a sense) . The resulting Ha — R dif- 
ference image, shown in Fig. is not perfect, partly be- 
cause of differences in the Ha to R flux ratios for different 
objects, and partly because of internal reflections in the 
Ha filter, which lead to low-level, out-of-focus images near 
every object. Nevertheless, the nebula can be seen more 
clearly, and in particular the confusion near the head is 
greatly reduced. The nebula has a cometary shape and is 
somewhat edge-brightened. 



3.1. Nebular spectrum 

In Fig. |l|, we show the spectrum of the nebula. 
This was formed by extracting a region of 2" around 
RX J1856.5— 3754 from the spectral images for both slits, 
subtracting the contribution from RX J1856.5— 3754, and 
averaging the two weighted by the exposure time. The 
flux calibration was done with the help of an exposure of 
the standard star EG 274 (Paper I). From the figure, it is 
clear that no other lines than Ha and H/3 are observed. 
The slightly elevated continuum at longer wavelengths re- 
flects an increased contribution of an extended source near 
the neutron star (see below). 



3.3. Orientation 

From Fig. |, it is clear that the nebula is aligned 
almost, but not precisely with the proper motion of 
RX J1856.5— 3754: the position angle of the nebular axis 
is 97° ± 1°, while the position angle of the proper motion 
is 100?3 ± 0?1 (Sect. 1). This shght but significant offset 
is not surprising, because the shape and orientation of the 
nebula is determined not by the motion of the neutron 
star relative to the Sun, but by its motion relative to the 
surrounding medium. 

For the ambient interstellar medium, it is reasonable to 
assume that its velocity is small relative to the local stan- 
dard of rest. If so, as observed from the Sun, the medium 
will have an apparent proper motion /i^P and radial veloc- 



4 



M. H. van Kcrkwijk and S. R. Kulkarni: An unusual Ha nebula around RX J1856.5— 3754 




« -56' 



ji- - 1 ' . 





^ •;***■' 

a;^! .'frifc^^ iP^I 





56^40= 



5Bm4os 



54-30" 



-54'45" 





56'na4s 





1 1 J 

^ ■ Ho-H 




* ' a 


1 1 


^ 1 . f T :, 



&6™36» 56™35s 5811345 ge^iaSii 56™35= 56™34» Se^ae^ &8™35» 56'"34' 5a'"36' 56'"35s 56™34s 

Fig. 2. Images of RX J1856.5— 3754, its nebula, and the field. The upper-left panel gives an overview of the stacked 
Ha image. In the upper-right panel, all but the most overexposed stars have been removed using point-spread function 
fitting. Except for the nebula (inside the box in the upper- left panel), all diffuse emission that is seen has matching 
emission in the B and R images, indicating that it is due to scattered light. The lower row of panels shows enlargements 
of the area inside the box in the upper-left panel, in B, R, Ha and Ha — R. In the Ha image, stars L, C, and F are 
indicated. For the Ha — R image, the R-band image has been optimally matched to and subtracted from the Ha image 
in order to remove the contribution from broad-band light. The tick marks indicate the position of RX J1856.5— 3754, 
as measured from the B-band image. In the R-band image, one sees that there are two objects close to it. One, about 
I'.'S to the ESE, has been detected as weU in the F606W WFPC2 images taken with HST ([Walter 200l|; star 115). 



The other, about 0'.'9 to the NE, is extended and only barely visible in the HST images. These objects give a false 
impression of a bright head of the nebula in the Ha image; this impression disappears in the Ha — R image. In the 
Ha — R image, the proper motion is indicated by the longer dashes, which connect the positions the object had between 
60 and 70 years ago and will have in 10 to 20 years from now. The dash marked '0' shows the proper motion as seen 
from the Sun, while that marked 'LSR' shows the proper motion as seen by an observer moving with the local standard 
of rest (and therewith moving with the interstellar medium local to the neutron star; see Sect. 3.3). 



ity "md. given by 



Mi' 



rad 



\U(7) sin / — 1^(7) cos I] 



ap 



- [{Uq cos I + Vq sin I) sin b - 



= — (C/q cos/ + T/q sin^) cos6 



- Wq cos b] , 
Wq sin b, 



(1) 



where (/,6) = (-1?40, -17?22) 
tude and latitude of RX J1856.5- 



are the Galactic longi- 

-3754, and (t/, V, W)q = 



(10.00, 5.25, 7.17) kms^i the velocit y of the Sun rela- 
tive to the local standard of rest (Dehnen & Binney 
1998). Inserting the numbers, we find 
(-19.3, 



-36.2)c?gQ masyr 
With the above. 



and v^P^ 



-7kms ^. 



the proper motion of 
RX J1856.5— 3754 relative to the nebula, and there- 
with, by assumption, relative to the local standard of 
rest, is (V^,Mb ) = (83.1, -289.7) masyr"! for d = 60pc. 
This implies a total proper motion fj}^'^ = 301.4 masyr^^ 
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at position angle 9 = 95? 3; the latter is close to the 
angle observed for the nebular axis (Fig. |). The implied 
relative spatial velocity on the sky is wjf — 86kms^^, 
while the relative radial velocity and total velocity, 
assuming an origin in the Upper Scorpius association, are 



,lsr 



,,lsr 



^ -50kms 1 and v'-^^ 



3.4. Flux calibration 



100 kms- 



No Ha calibration images of emission-line objects were 
taken, and therefore we used our carefully calibrated spec- 
tra of stars F, L, and C (Paper I) to estimate the through- 
put. For this purpose, we measured instru mental magni - 
tudes on our Ha images, using daophot ( Stetson 1987 ) 
in the same way as was done for B and R in Paper I, and 
compared these with the expected magnitudes, 

= -2.51ogy' AuTi?a(A)i?opc„(A)-^/A(A)dA, (2) 4. A pulsar bow shock? 



6628 A. The H/3 and H7 surface brightnesses were deter- 
mined similarly. Note that the Ha surface brightnesses 
inferred from the spectra are more certain than those in- 
ferred from the Ha image, since they suffer much less from 
uncertainties in the background level. 

The nebula is clearly resolved in all directions. It 
reaches half-maximum brightness I'.'O ± 0'.'2 in front of the 
neutron star, about 2" to each side, and ~ 3" behind the 
neutron star. The tail is visible to at least 25" behind the 
neutron star. The peak surface brightness of the nebula 
is 1.3 X lO"'^^ ergs"^ cm~^ arcsec"^, corresponding to a 
photon rate of 4.3 x 10^^ s~^ cm~^ arcsec"^. 

Below we discuss the origin of this nebula and see if 
we can use it to shed additional light on the nature of the 
enigmatic object, RX J1856.5— 3754. 



where Aut = 51.2 m^ is the telescope area, Ra is the re- 
sponse of the Ha filter, i?opon the combined response of all 
other elements (atmosphere, telescope, collimator, detec- 
tor), and /a (A) the spectrum of the object in wavelength 
units. For Ra, we use the response curve from the ESO 
web site|J] as measured in the instrument using a grism; 
the filter has a peak throughput of ^ 65% and full width 
at half maximum of 60 A. We infer a total efficiency for 
all other components of i?open = 32 ± 1% (assuming i?opcii 
does not vary strongly over the Ha bandpass; we verified 
that a similar analysis for the R band gives a consistent 
result). Here, the uncertainty includes in quadrature a 2% 
uncertainty in the flux calibration of stars L, C, and F. 
For the conversion from Ha count rates to photon rates, 
therefore, only the filter efficiency at Ha is required. We 
find Ra{Xa) = 60 ± 5%, where the relatively large un- 
certainty stems from the fact that there is some depen- 
dence of the filter characteristics on position on the detec- 
tor (even though the interference filters are located in the 
convergent beam). 



Cometary nebulae shining in the Balmer lines (with no 
detectable lines from any other element) have been seen 
around a number of pulsars. The accepted physical model 
for such nebulae is the pulsar bow-shock mo del, which 
was invoked first by Kulkarni fc Hester (1988 ) to explain 
a similar Ha nebula around the binary millisecond pulsar 
PSR B1957+20. Thus, both the similarity of the nebula to 
nebulae seen around some pulsars, as well as the reason- 
able possibility that RX J1856.5— 3754 is a million-year 
old neutron star (see Sect. ^) and thus might be a pulsar, 
motivate us to consider the bow-shock model. 

4.1. Stand-off distance and shape 

Consider a pulsar moving supersonically through partially 
ionised interstellar medium of density pism- Along the di- 
rection of motion the medium ram pressure is pism^ng, 
where v^s is the speed of the pulsar. At the apex of the 
nebula, the (quasi) spherical relativistic pulsar wind is bal- 
anced by this ram pressure at a "stand-off" distance ra, 
where 



3.5. Results 

In Fig. m we show the calibrated Ha — R image, rotated 
such that its axis of symmetry is horizontal (assuming a 
position angle of 97°, see Sect. ^.3[ ). Note that the level of 
the background remains somewhat uncertain; indeed, it is 
the limiting factor in the resulting Ha fluxes. Also shown 
in the Figure are cross-cuts through the image, one along 
the axis of symmetry, and four across it at different dis- 
tances from the head of the nebula. Furthermore, the Ha, 
H/3, and H7 surface brightnesses inferred from the spec- 
tra are shown. The Ha brightnesses were derived from the 
spectral images, by integrating over the wavelength range 
of 6554-6570 A and subtracting the average of the inte- 
grated fluxes in wavelength ranges 6496-6550 and 6574- 



E 



^ |http : //www, eso . org/ instruments/for si /Filters/ 
FORS_filter_curves.html, filter halpha+59. 



A'KrIc 



(3) 



Here, E is the luminosity in the relativistic wind, which is 
assumed to be primarily in the form of relativistic particles 
and magnetic fields. B y considering momentum balance 
along other directions, Wilkin (1996| ) derived an analytic 
solution for the shape of the nebula. 



Tg sin 9 




tan 9 



(4) 



Note that in this derivation it is assumed that the shocked 
gas cools rapidly, so that thermal pressure is negligible (a 
two-dimensional equivalent to the 'snow-plough' phase in 
supernova remnants). If this assumption is incorrect then 
the tail should be slightly wider. 
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Fig. 3. Flux-calibrated Ha images and cross cuts of the nebula around RX J1856.5— 3754. The lower-left panel shows 
the nebula, using coordinates relative to the neutron star, vifith the horizontal axis representing the major axis of 
symmetry of the nebula. In the upper-left panel, the surface brightnesses along the axis of symmetry of the nebula are 
shown. The points with error bars represent the Ha, H/3, and H7 brightnesses inferred from the spectra taken with the 
slit positioned over star F, very close to the axis of symmetry; these are marked as 'a,' '/?,' and '7,' respectively. See 
Fig. 2 and Fig. 1 in Paper I for star designations including star F. The nebula is faint in H/3 and undetectable in H7. 
All points are averages over 3-pixel (0'.'6) intervals along the slit. The thin lines indicate the zero level and show the 
point-spread function of star F (also averaged over 3 pixels). For Ha, it also shown at the position of RX J1856.5— 3754. 
The histogramme represents Ha measurements inferred from the Ha — R image, by averaging over the numerical slit 
marked 'z' in the upper-right panel; its zero level is again indicated by a thin line. The grey scale used in the images 
is next to the upper-left and above the lower-right panels. In the lower-right panel, the cross cuts perpendicular to the 
axis of symmetry are shown. The points with error bars represent the brightnesses inferred from the spectra taken with 
the slit positioned over star L, and the histogrammcs those inferred from the image, by integrating along numerical 
slits marked 'p,' 'q,' 'r,' and 's.' The zero levels and point-spread functions arc indicated by thin lines. Also displayed 
in the upper-right panel is the expected shape for a bow-shock (see Sect. H). 



4.2. Shock emission mechanisms 

Along the cometary surface, one expects two shocks, with 
in between an inner region consisting of shocked pulsar 
wind and an outer region of shocked interstellar medium. 
The former will shine primarily via synchrotron emission 



and is best searched for at X-ray (e.g. the Mouse, Predehl 



fe Kulkarni 1995) or radio wavelengths (see Gaenslcr et al, 
200c for a review of radio searches). 



All astrophysical shocks appear to be coUisionless and 
mediated via ionised particles. Thus, the ambient neutral 
particles (i.e., atoms) do not notice the shock, but do find 
themselves subjected to collisional excitation and ionisa- 
tion by shocked electrons and to charge exchange with 
shock ed protons ( phevalier fc Raymond 1978| ; Raymond 
1991 ). Both these processes can lead to line emission 



Electronic excitation of the unshockcd neutrals will result 



ing the velocity of the shocked proton and its Ha emission 
will reflect the large velocity dispersion and the bulk ve- 



locity of the shocked protons (Chevalier & Raymond 1978 



Raymond 1991). Elements other than Hydrogen will also 
be excited coUisionally, but the large abundance of hy- 
drogen and the high temperature ensure that the emitted 
optical spectrum is dominated by hydrogen lines. Hence, 
the name "Balmer-dominated" shocks. 

Balmer-dominated shocks are rare because two strin- 
gent conditions need to be met. First, the ambient ISM 
has to be (partially) neutral. Second, the shock has to be 
"young." The latter condition arises from the fact that for 
shock speeds greater than 110 kms""'^, the post-shock gas, 
if allowed suflicient time to cool, produces so much ion- 
ising flux that the pre-shocked ambient gas ahea d of the 
shock front becomes almost completely ionised ( ^huU fc 



in a n 4rrow line, reflecting the velocity dispersion of the McKee 1979 ; McKee fc HoUcnbach 198C ). Thus, B aimer 



unshocked gas. Charge exchange of a Hydrogen atom with 
a fast-moving proton will result in a Hydrogen atom hav- 



dominated shock spectra are expected only in shocks 
which have not enough time to sweep a "cooling column" 
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density, A^cooi = l.SxlO^^vfi^l cm-^ (Hartiganct al. 1987) 



here, wioo is the shock speed in units of 100 km s ^ 

Given this reasoning, it might come as a surprise that 
Balmer-dominated bow shocks have been found around 



pulsars with a wide variety of ages (Kulkarni & Hester 



1988; (hordes et al. 1993; BeU et al. 1995). Pulsar bow 



shocks are far from planar, however, and as a result the 
shocked material subtends an ever decreasing solid angle 
as seen from the vantage point of the stand-off region. In 
addition, pulsars spend most of their time moving through 
diffuse interstellar medium and thus matter recombines 
only far behind the pulsar. The combination of both rea- 
sons ensures Balmer-dominated pulsar bow shocks can 
form. 



4.3. Application to RX J1856. 5-3754 

Turning now to RX J1856.5— 3754, we start by noting that 
the shape of the nebula fits reasonably well with the ex- 
pression given by Eq. ^, as can be seen in Fig. By try- 
ing different inclinations i and projected angular stand-off 
distances , we find a best match around i — 60° and 
^11 = I'.'O, with a 'chi-by-eye' uncertainty of about 15° and 
0'.'2, respectively. This determination is consistent with 
our inference of the radial velocity with respect to the 
local medium (see Sect. H). Note, however, that, as men- 
tioned above, Eq. ^ is derived under the assumption of 
rapid radiative cooling; if this is not met, the tail would 
be wider and hence the observations would imply inclina- 
tions closer to 90°. We return to this aspect below. 

In the framework of the bow-shock model, apart from 
the orientation, there are four unknowns, the relative 
velocity, the pulsar wind energy loss E, the interstellar 
medium density pism, and the fraction = nno/nn of the 
neutral ambient gas. In contrast, we have only three mea- 
surements: the velocity, the observed Ha flux, Ua, and the 
angular distance between the neutron star and the apex of 
the cometary nebula, . The physical length of the stand- 
off distance is ra = d0\\/ sini = GOdeo/ sinz AU, where dgo 
is the distance in units of 60 pc. For simplicity, we will 
drop the dependence on sinj in our estimates below, and 
will take Vns = d/ias = lOOdeo kms""'^. 

The Ha flux is proportional to the incoming flux of 
neutral particles (H i atoms) 



Raymond (1991) states that 



for a broad range of shock velocities, about 0.2 Ha pho- 
tons are produced per neutral particle before that particle 
is ionised in the post-shock region. Unfortunately, there 
are few calculation for the yield of Ha, the number of 
Ha photons per incoming neutral atom, fa, for low veloc 
ity shocks, ~ 100 km s 
[(1985| ), for a 120 km s 



. According to Cox 
shock, 



Raymond 



fa — 0.1 for case A (i.e., 
most higher order Lyman photons, Ly/3, Ly7, etc., escape 
from the shocked region) or ~ 0.23 for case B (higher order 
Lyman photons are absorbed locally and re-emitted until 
they are converted; e.g., to Ha and two-photon continuum 
for Ly/3). We will assume case A and normalise the yield 
accordingly: fa = 0.1/a^o.i- 



Considering the head of the nebula only, the expected 
number of Ha photons is related to the other parameters 
as follows: 



/ariHOTrr^Wns _ 1 



47rd2 



-faniio0^^1nsd. 



(5) 



Here, r± is the radius of the bow shock perpendicular to 
the axis of symmetry and 6± = r±/d the corresponding 
angular distance. 

With 9± = 2", the measured /ins (Sect. H) and the 
observed Ha photon rate of Ua 
find 



rT-HO = 0.8 cm ^ fal.idw- 



(6) 



Noting that the total density of the ambient gas is 
Pism = '^howh-^~^^o'^: where X = 0.73 is the fraction of 
hydrogen by mass and = 'in" /"-h is the neutral fraction, 
and applying Eq. ^, we find. 



i?~6xl03iergs-i/alieo"'di 



-U3 
60- 



(7) 



We find that because we have fewer measurements than 
unknowns, we are only able to constrain the quantity ^o^'- 
Since Co ^ Ij we obtain a lower limit to E. 

If, as argued above, RX J1856.5-3754 is a t = lO^g yr 
old pulsar, then assuming that the main energy loss is due 
to magnetic dipole radiation (magnetic field strength B, 
braking index n = 3), we can apply the usual relations, 
Tchar = P/2P = t/x, with X = 1 - [Po / Pf - 1 (whcrc 
P and Pq are the current and birth period, respectively), 
E = Mtt^P/P^ = J47rV2Tp2 (where / = 10''^/45gcm2 
is the moment of inertia), and B = 3.2 x 10^^{PP)^^'^ = 
3.2 X 101V/(2t)1/2 G, to obtain. 



/2 ,3/2 



P ~ S 3s y1/2/1/2,-1/2 .1/2 .1/2,: 

B ^ l.5xlO'^GxlU\'f'JLeo^^d] 



'60 

1/2^3/2 
60 • 



(8) 
(9) 



Since x ^ ^ and < 1, the values of B and P are true up- 
per limits. The limits on B and P are consistent with those 
expected from ordinary million-year old pulsars. Thus, the 
bow-shock model appears to provide a reasonable descrip- 
tion of the nebula we observe. 

4.4. Consistency checks 

We now verify that the bow-shock model is consistent with 
the assumptions we made. First, the shape of the bow 
shock was assumed to be specified by Eq. || and this as- 
sumption implicitly requires that the shock radiate much 
of the energy. How good is this assumption? Every Hi 
atom comes into the shocked region with an energy of 
^mni^ns — ^^"^ns 100 where Wns.ioo is the shock speed 
in units of 100 kms^^. Of this, about jq goes into shock 
heating and the remainder into bulk motion. Only one out 
of five hydrogen atom are excited to the n = 3 level before 
being ionised; including other levels, likely the energy loss 
from the shocked region is ~2 x 13.6 eV per neutral. Thus 
the ratio of energy radiated to the incoming kinetic energy 
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is ^ l-O^o- Thus, the use of a "snow-plough" approxima- 
tion is not unreasonable. This is supported observation- 
ally by the fact that the three pulsar bow shocks studied 
to date, in particular the two with velocities comparable 
to tha t of RX J 1856. 5-3754 (PSR B1957+20, [Kulkarn 
& Hes fer (498^ ); PSR J0437-4715, [Bell et al. (I995| )) 
display nebulae with sharp boundaries. 

Second, the ambient gas has to be partially neutral in 
order to see Balmer line emission. The inferred Hydrogen 
number density of ~ f cm~'^ is not an unreasonable den- 
sity (in order) for the Warm Neutral Medium (WNM) 
or the Warm Ionised Medium (WIM). The two phases 
taken together are expected to occupy 50% of interstellar 
space locally; see Kulkarni fc Heiles (4988D and Dickey & 
Lockm m (4990D for general reviews of the conditions in 
the interstellar medium. The ionisation in the WNM is low 
but th e W I M is partially ionised, 
et al. l |995t [Rcdfield fc Linsky 2000]) . 



0.5 (see Reynolds 



Independent of the nature of the ambient gas, the 
bow-shock model requires that the incoming gas not be 
fully ionised by the extreme ultra-violet radiation from 
RX J1856.5— 3754. In the next section, we treat the case 
of a pure ionisation nebula. We find that the ambient gas 
is becomes half- ionised at an angular offset, in the for- 
ward direction of only '-^O'.'2/bb (Eq. |2^) where /bb is the 
excess of ionising flux (in the extreme ultraviolet band) 
over that expected from blackbody fits to the X-ray and 
optical data of RX J1856.5— 3754. Thus we would have to 
abandon the pulsar bow shock model if /bb ^ 4. 

Third, in our estimates, we implicitly assumed the 
bow-shock was thin compared to the distance to the 
neutron star. Ghavamian et al. (2001) have computed 



the ionisation structure of Balmer-dominated shocks for 
riHO = lcm~'^ and = 0.5. For a 250kms~^ shock, 
they find that the time scale on which an H i is ionised by 
electrons is 0.2 yr, and that the time scale on which pro- 
ton charge exchange becomes ineffective is about ~0.4yr. 
These time scales are expected to depend only weakly on 
velocity and to be inversely proportional to density. The 
corresponding length scales are I ~ 10"'^'*Wns,ioo cm, which 
is smaller than ~ lO^^dgo cm. Thus, we conclude that 
the assumption of a thin bow shock is reasonable. 

Finally, we verify whether the assumption of case A 
was reasonable. For the Lyman lines, the line-c entre cross- 
section is given by ( Rybicki fc Lightman 1979| ) 



continuum), we conclude that our assumption of case A is 
reasonable. 



5. An ionisation nebula? 

RX J1856.5— 3754 is a source of ionising photons and, 
as with hot white dwarfs, should have an ionised region 
around it. Such a region would produce Ha emission, and 
would be deformed b ecause of the prop er motion of the 
neutron star. Indeed, Blaes et al. (1995| ) showed that iso- 
lated neutron stars accreting from the interstellar medium 
should have associated Ha nebulae with cometary shape. 

To see whether ionisation could be important for the 
nebula around RX J1856.5— 3754, we will first make an- 
alytical estimates in Sect. 



5.1, following the analysis of 



Blaes et al. (1995). We will find that ionisation could 



lead to a nebula with properties in qualitative agreement 
with the observations. The inferred conditions in the neb- 
ula, however, are rather different from those encountered 
in "normal" ionisation nebulae; we discuss this briefly in 
Next, in Sect. 5.3, we present a detailed simula- 



Sect. 5.2 



tion with which we attempt to reproduce the observations 
quantitatively. 

5.1. Analytical estimates 

In general, the properties of ionisation nebulae are deter- 
mined by the balance between ionisation, recombination, 
heating and cooling. In order to determine the relative im- 
portance of these processes, we first compare their charac- 
teristic timescales with the crossing time, i.e., the time in 
which the neutron star with velocity v^s moves a distance 
equal to the distance fa to the apex of the nebula, 



f — —— — —— — 3vr 

''Cross ^ y ■ 

^ns f^ns 



(11) 



For our analytical estimates here and below we take ra = 
d^a — 9 X lO^^dgQ cm and = rf^ns — lOOdgo kms~^, i.e., 
we ignore the effects of changes in inclination and the dif- 
ference between observed proper motion and motion rela- 
tive to the local interstellar medium discussed in Sect. |3.3[ 
We will also assume pure Hydrogen gas. We will include 
Helium and use the correct geometry in the numerical 
model described in Sect. 5.2. 



— f / "^H 

rricC uq V 27rfei 

= 1.16X 10-"cni2 XoAh,nT^'^\ 



(10) 



where /i.„ is the oscillator strength (for Ly/3 through 
Lye, /i,2 5 = 0.4162, 0.07910, 0.002899, and 0.001394). 
The line-centre optical depth in the Ly/3 line is thus 
~2 X a^g^i^iuiiol ~ 0.8; here the multiplicative factor of 2 
is a crude approximation to account for the bow shock ge- 
ometry. Since in every absorption, there is a probability of 
only 0.118 that Ly(3 is converted to Ha (plus two-photon 



5.1.1. Ionisation and recombination timescales 

The typical time for a neutral particle at the apex of the 
bow shock to be ionised by the emission from the neutron 
star is 



t-ioii ~^ \ A n ^ii 

\ATrrl 



(12) 



where A^x is the rate at which ionising photons are emit- 
ted, nx = Nx/4:TT(P the rate at which these would arrive 
at Earth in the absence of interstellar absorption, and ajon 
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the effective cross section for ionisation of Hydrogen. The 
latter is given by 



1.4 X 10-19 cm^T-ig^T. (^3) 



here 5*,^ is the surface emission, which for the numerical 
estimate is assumed to be a black-body with temperature 
fcToff = 50Tns,5o eV. For this cross section, the effective 
optical depth r = nnoaion'^a — 1-3 x lO-^nno^ is very 
small (here njjo ~ 1"-ho i cm"^ is the neutral Hydrogen 
number density). Even the optical depth at the Lyman 
edge is small: acdgc = 6.3 x IQ-^* cm^ and hence Todgc — 
6 X IQ-^nHO,!. 

We estimate nx by scaling to the observed, de- 
reddened optical flux, and assuming that the emitted spec- 
trum resembles that of a black body, i.e.. 



^opt 



0.4s" 



/bbTj 



ns,50' 



(14) 



where we used the unabsorbed photon rate inferred in 
Paper I from the optical and ultra-violet photometry of 
8.5 X lO-^s-^cm-^ at 5000 A and where /bb is a 
factor which takes into account the extent to which the 
emitted spectrum deviates from that of a black body. 
With the above estimates, we find 



Uon ~ 13yr fy^^^T^ 



0.33 
bb -^ns,50 • 



(15) 



The typical time it takes for a proton to recombine is 



iroc = (arccJ^c) ^ ^ 10^ yr "o 1 ' 



(16) 



where a^cc is the recombination cross section and no = 
Irie,! cm"^ the electron number density. For the estimate, 
we used aroc — 4 10^"'^'^ cm^^ s^^, which is for case A at 
lO** K (Osterbrock 1989). Case A is appropriate here, since 



we found above that the Lyman continuum is optically 
thin and thus recombination directly to the ground state 
does not lead to local ionisation of another atom (as is 
assumed in Case B, for which the recombination rate is 
~40% lower). We will find below that the temperature is 
likely substantially higher than 10* K, but this will only 
lead to lower recombination rates. 

Clearly, for any reasonable density, the recombination 
time scale is far longer than the crossing time. The ion- 
isation time scale, however, is comparable to the cross- 
ing time. To match the two, /bb will have to be some- 
what larger than unity, but such a deviation may simply 
reflect the difference between a real neutron star atmo- 
sphere and a black body spectrum. Thus, the size of the 
nebula is consistent with it being due to ionisation. The 
tail of ionised matter will have length of /ins^rcc — 10°. 
This is far longe r than the observed length, but we will 



see in Sect. 5.1.3 , where we consider thermal balance and 
the relevant emission processes, that this difference can 
be understood. First, however, we calculate the ionisation 
structure of the head of the ionisation nebula. 



5.1.2. Structure of the ionisation nebula 

In general, the fraction of ionised particles ^+ = n^/nn 
changes as a function of time as 



dt 



Nxe- 



-T(r) 



47^7-2 



(17) 



Given our estimates above, both optical depth effects and 
recombination can be ignored. For this case, the ionisation 
structure can be solved analytically. In a cylindrical coor- 
dinate system (p, z) with the neutron star at the origin 
and gas moving by in the +z direction, the solution is 



^+ = 1 — exp arctan — 

p —z 



where the length scale zq is given by 



zq 



47rw„s 



(18) 



(19) 



and where arctan(p/— z) - taken to be in the range [0,7r] 
- measures the angle with respect to the direction of the 
proper motion. From Eq. one sees that in a parcel of 
interstellar matter at given impact parameter p, ^+ in- 
creases slowly at first, then more rapidly as it passes the 
neutron star, and finally slowly approaches its asymptotic 
value of 1 — e~'^^°/P. For large p, ^+ oc 1/p, as expected 
when integrating over a 1/r^ distribution. 

The angular size corresponding to the length scale zq 

is 



d 



eS 
f^ns 



• ^ /bbJns,50 



(20) 



Like above, we find that to match the observed size, a 
somewhat larger ionising photon rate appears required 
than expected based on the observed optical and X-ray 
fiuxes, i.e., /bb > 1- We will return to this in Sect. |3.3| , 
where we describe our more detailed modelling. 

5.1.3. Thermal balance and emission process 

We first ignore all cooling and calculate the temperature 
due to heating by the photo-ionisation. Using a typical 
energy e of the ionising photons of 



hvau{Sv/hv) dv 



J^^^ ay[Su/hv) dv 
we find for a pure Hydrogen gas 



31.2 evr; 



0.25 
ns,50 ' 



(21) 



where Xi.oo = 13.6 eV is the ionisation potential of hydro- 
gen. 

The high inferred temperature of ~5 x lO** K for ^_|_ = 
0.5 has an important consequence, viz., that coUisional ex- 
citation and ionisation of Hydrogen are important. These 
processes will dominate the cooling. Furthermore, the Ha 
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emission will be dominated by radiative decay of collision- 
ally excited neutrals. Below, we first estimate the density 
required to produce the observed Ha brightness. Next, we 
use this to estimate the cooling time scale and check con- 
sistency. 

The rate for coUisional excitation to the n = 3 level 
is nenHf)'7i,3 = ^+ (1 — i+)'n'^qi,3, w ith the rate coefficient 



9i,3 given by (see Ostcrbrock 1989) 

e-xi,3/fcT 



91,3 — ^1,3 



2y/TracaQ 



Xl,c 



kT 



(23) 



where Ti,3 is the sum of collision strengths to the 3s, 
3p, and 3d states, uji — 2 the statistical weight of the 
ground state, xi.3 — 12.1 eV the excitation potential, and 
20raca^ = 2.1716 x lO^'^cm^s"^ The collision strength 
varies slowly, from T1.3 = 0.25 at 10^ K to 0.43 at 7 



10* K ( Anderson et al. 2000 ). Inserting numbers, one finds 



that the largest volume emission rates are produced where 
^+ ~ 0.7. In consequence, one expects that the emission 
will peak at some distance from the neutron star, and 
that the ionisation nebula will have a somewhat hollow 
appearance. 

In order to estimate the hydrogen number density re- 
quired to reproduce the observed brightness, we will as- 
sume case A, i.e., the medium is optically thin to Ly/? 
photons, so that of the excitations to the 3p state only 
11.8% lead to emission of an Ha photon (for 3s and 3d, 
the transition to the ground state is forbidden, so all exci- 
tations lead to Ha emission). We label the appropriately 
corrected excitation rate as q^. Furthermore, we will ig- 
nore excitations into higher levels which lead to Ha pho- 
tons. We consider the total Ha photon rate integrated 
over the slit (with width Wsia — 1") perpendicular to the 
proper motion. The expected rate is given by 



7rrj^dwsiit»H""'cga 
47rd2 

-9']_dwsiit£.+ {l - Mnfiq, 
1.2 X 10 



2 „A 



(24) 



where for the numerical estimate we used ^+ — 0.5 and 
q^ ~ 2.2 X 10~i° cm3 s^^ (appropriate for T = 5 x 10"* K). 
To match the observed photon rate of 1.7 x 10^^ s^^ cm"^ 
along the slit, thus requires 



4 cm ■^d, 



-1/2 



60 



(25) 



This is somewhat denser than the typical density of the in- 
terstellar medium or the mean density of N-n/d ~ 1 cm~^ 
along the line of sight, but not by a large factor. Hence, the 
observed brightness can be reproduced without recourse 
to extreme assumptions about the interstellar medium in 
which the source moves. 

With the above estimate of the density, we can cal- 
culate the cooling time scale and verify a posteriori that 
cooling is not important near the head of the nebula. As 
mentioned, at the high inferred temperatures, the domi- 
nant cooling processes are collisional excitation and ionisa- 
tion of Hydrogen. At T ~ 5 x 10* K, collisional excitation 



of the n = 2 level is most important, with a rate coef- 
ficient (71,2 ~ 2 X 10~^cm=^s~^ (Eq. |2|, xi.2 = 10.2 eV, 



T1.2 = 1.4; Anderson et al. 2000). Hence, the cooling time 
scale tcooi is roughly 



cool 



Xi,. 



(nn + nc)\kT _ 

rtH"'T-cgi,2Xl,2 (1 - ■?+)?^H'?l,2Xl,2 

14 yr n^\. 



(26) 

This time scale is a little longer than the crossing time, 
which confirms that temperatures as high as those inferred 
above can be reached. The time scale is sufficiently short, 
however, that more precise estimates will need to take 
cooling into account. We will do this in Sect. |5^ . 

The relatively short cooling time scale provides the 
promised answer to the puzzle posed in Sect. ^.l.l| : the tail 
appears far shorter than the recombination length of 10° 
because the temperature and hence the emissivity drop 
rapidly behind the neutron star. In principle, the tail will 
emit recombination radiation, but this will be at unob- 
servable levels. 

With the inferred temperatures, we can also under- 
stand why we see Balmer emission only (Fig. For 
Helium, the temperature is too low for collisional exci- 
tation to be important. This leaves recombination and 
continuum pumping, neither of which are efficient enough 
to pr oduce observable emission in our spectrum (see also 
Sect. ^.3.4 ). For the metals, collisional excitation may well 
take place, with rate coefficients comparable to or even 
larger than those for Hydrogen, but given their low abun- 
dances relative to Hydrogen, the resulting emission will be 
unobservable as well. 



5.1.4. Dynamical effects 

The heating due to photo-ionisation by the relatively hard 
photons from the neutron star will also lead to an over- 
pressure in the gas near the neutron star relative to that 
further away. As a result, the gas should expand. We es- 
timate the velocity in the direction perpendicular to the 
proper motion from Vp = Opt, where Op is the accelera- 
tion associated with the pressure gradient perpendicular 
to the proper motion and t is the time scale on which the 
pressure gradient changes. The acceleration is given by 

IdP 2 d^^ 

Up — — = ( f — Vi -.1 

P ap 



3mH 



(.-Xl,oo)^ 



(27) 



At z = 0, one has d^+/dp = —{po/p'^)e~P°^P, where po = 
■|zo- The pressure gradient changes due to expansion and 
due to the neutron star passing by. The time scale for the 
latter is much faster, and is approximately given by 

~Po /p 

.. (28) 

z=0 



1 

t ~ 



d^+/dz 



Thus, for the velocity one finds 
2(e-Xi,oo)po 



llkms-i 



T, 



0.45 



ns,50 "60 ' 



(29) 
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where the numerical value is for p = pq. Note that for large 
p, the expansion velocity becomes independent of p. This 
is confirmed by our detailed models, but only if cooling is 
unimportant. 

Because of the expansion, behind the neutron star the 
gas will be more tenuous and emit fewer Ha photons. This, 
in turn, leads to a more hollow appearance of the emis- 
sion nebula as a whole. We include the expansion in our 
numerical model, which we describe in Sect. 5.2. First, 



however, we discuss briefly why the properties we derive 
for the nebula differ greatly from those of "normal" ioni- 
sation nebulae. 



5.2. Intermezzo: comparison with "normal" ionisation 
nebulae 

Ionisation nebulae observed around hot stars, such as 
planetary nebulae and H ii regions, typically show a wealth 
of emission lines of many elements, and have inferred tem- 
peratures around 10^ K ( psterbrock 1989| ). In these, the 
medium is in (rough) equilibrium: ionisation is balanced 
by recombination and heating by cooling. The emission 
from Hydrogen is due to (inefficient) recombination, while 
that of the metals is due to (efficient) collisional excitation; 
hence, the metal lines are strong despite the low metal 
abundances. The emission is strongest where the matter 
is densest, and nebulae with roughly uniform density ap- 
pear filled (e.g., Hii regions). 

The situation for RX J1856.5— 3754 is completely dif- 
ferent: there is no equilibrium; recombination is not rele- 
vant at all and cooling only marginally so. The reason is 
that the neutron star emits only few ionising photons (be- 
cause it is small) and moves fast. As a result, the nebula 
is "matter bounded," and much smaller than would be in- 
ferred from the Stromgren, "ionisation bounded" approx- 
imatiorj^ (see Blaes ct al. 1995). The interstellar medium 
simply does not notice the few photons passing by until 
the neutron star is very close, at which time the photo- 
ionisation and heating rates increase far more rapidly with 
time than can be balanced by recombination and cooling. 

In some sense, in considering the conditions in the ioni- 
sation nebula derived here, it may be more fruitful to com- 
pare with the conditions in a bow shock rather than with 
those in "normal" ionisation nebulae. Like for the nebula 
considered here, in the case of a bow shock, the ambi- 
ent medium receives a sudden injection of energy, which 
cannot be balanced by cooling, and results in a very high 
temperature (Sect. ^). In both cases, collisional excitation 
and ionisation of Hydrogen are highly important, leading 
to Balmer-dominated emission spectra. The main differ- 
ence is that for the case of a bow shock, not only energy 
but also momentum is injected into the medium. The ab- 
sence of the shock associated with the momentum injec- 



^ The i ncorrect applicatio n of the Stromgren approximation 
likely led Pons et al. (2001 ) to their remark that the nebula 



tion makes the case of the ionisation nebula a much more 
tractable one. 



5.3. A detailed model 

Our analytic estimates showed that the nebula could be 
due to ionisation, without requiring extreme conditions. 
To make a more quantitative comparison, we made a 



model along the lines of that described by Blaes et al 
|(1995| ), except that we include the dynamics. For this pur 



pose, we use the hydrodynamics code zeus-2d (Stone & 
[Norman 1992 ). We changed the code to allow it to keep 
track of abundances of separate ions, and wrote a routine 
that calculates the photo- ionisation rates, the associated 
heating, and the cooling due to collisional processes. We 
only include Hydrogen and Helium (at cosmic abundances, 
X = 0.73 and Y = 0.27); metals are not important for the 
ionisation balance, and, as argued above, do not have their 
usual role as coolant in our case. 

In our calculations, we take into account ionisation 
of Hi, He I, and Hen, as well as the associated heating. 
For each of the ions, we calculate effective cross sections 
and energies e as in Eqs. 13 and Ell by integrating 



a black body photon spectrum over the appropriate cross 



sections (Hi and Heii: Osterbrock 1989; Hei: Verner et al 



1996). We verified that optical depth effects were negligi- 



ble. Note that for hard ionising spectra, Helium is ionised 
before Hydrogen, and that the resulting photo-electrons 
have larger energy (for kTns = 50 eV: e — Xi.oo — 18 eV 
for Hydrogen and ~ 33 eV for both He i and He ii) . As a 
result, temperatures are high already when the fraction of 
neutral Hydrogen is still relatively high, and the nebula 
will be brighter at large distances than would be expected 
for pure Hydrogen. 

For the cooling, we include collisional ionisation of 



Hydrogen (Scholz & Walters 1991) as well as collisional 
excitation of Hydrogen from the ground state to levels 2- 



5 (Anderson ct al. 200C). We verified that recombination 
was unimportant; free- free radiation is also negligible. We 
considered secondary ionisations by the fast ionised elec- 
trons, using the values found by 3hull (1979 ), but found 
that it was unimportant in any region with ^+ > 0.05 and 
hence had only very slight overall effect. For the sake of 
simplicity, therefore, we have not included it in the results 
presented here. 

In order to compare our results with the observations, 
we use the results of the simulation to calculate the colli- 
sional excitation rates from the ground state to all in divid- 
ual levels 2s to 5g of Hydrogen ( Anderson et al. 200C ), and 
then infer emission in particular transitions using the cas- 
cade matrix. As before, we assume that all Lyman lin es are 
optically thin (case A); we will return to this in Sect. 5.3.4 
below. 



around RX J1856.5— 3754 could not be due to ionisation. 
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Fig. 4. Model for an ionisation nebula. The coordinate frame is cylindrical, with z running along the axis of symmetry 
and p perpendicular to it. The neutron star is at the origin and moves to the left. The computational domain was 
between —18 x 10^^ < z < 36 x 10^^ cm and < p < 18 x 10^^ cm. The three left-hand panels show predicted runs 
of the fractional abundances of Hi, Hei, and Heii. Darker shades indicate lower fractional abundance; the 20, 40, 60, 
and 80% levels are indicated by contours. The right three panels show runs of density, temperature, and Ha emissivity. 
For the density, darker shade indicates lower density; the white contours are drawn at 3, 4, 5, 6, and 7 x lO"^"* g cm~'^. 
Overlaid are black contours and arrows indicating the velocity; the contours are drawn at 2, 4, 6, and 8kms~^. For 
the temperature, darker shade indicates higher temperature; the black contours are drawn at temperature of 1, 1.5, 
2, and 2.5 x 10^ K, and the white ones at 3, 4, and 5 x 10* K. Finally, for the Ha emissivity, darker shade indicates 
stronger emissivity; the contours levels are at 1/8, 1/4, 1/2, 1 (black), 2, 4, and 8 x lO^^*' s^^ cm^"^ (white). 



5.3.1. Model parameters 

Our model requires ten parameters, which describe the 
conditions in the undisturbed interstellar medium (1-5; 
density, temperature, and fractional Hi, Hei, and Hen 
abundances), the ionising radiation field (6, 7: neutron- 
star temperature and radius), the relative velocity be- 
tween the neutron star and the medium (8), and the geom- 
etry relative to the observer (9, 10: distance and inclina- 
tion). All but one of these, however, are constrained by ob- 
servations. The distance is around 60 pc, and the inclina- 
tion about 60° assuming an origin in Upper Sco (Sect, p); 
the implied relative velocity is lOOkms"^ (Sect. 3.3). 
Furthermore, the temperature of the neutron star is ^ 
50eV and the effective radius 7.0kmd6o(/bb/Tiis,5o)^^^j 
where /bb measures the deviation of the spectrum from 
that of a black body (Eq. |l^) and should be close to unity. 
Finally, the interstellar medium can be in one of four 
phases, the so-called cold-neutral, warm-neutral, warm- 



ionised, and hot-ionised phases (see Kulkarni & Heilcs 
IQSSfl^r a review). Of these, the hot ionised phases can- 
not be applicable to our case, since these would not lead 
to any emission, while the cold-neutral phase is unlikely 
because of its small filling factor. Furthermore, the warm 
ionised medium, with ^+ ~ 0.5 (Sect. H) is unlikely, since 



it typically is less dense than required. Thus, we will as- 
sume the medium is in the warm- neutral phase, i.e., it is 
neutral and has T ~ 8 x lO'^K. It should have a particle 
density in the range 0.1 to 10 cm^''. 

Given the above constraints, the only unknowns are 
the interstellar medium density pism and the ionising pho- 
ton rate factor /bb- As can be seen from the analytical 
estimates, these two have rather different effects: the den- 
sity sets the brightness of the nebula, while the ionising 
photon rate sets the scale. 

5.3.2. Model results 

In Fig. U, we show the results of a model calculation with 
Pism and /bb chosen to roughly reproduce the observa- 
tions: /bb — 1-7 and pism = 7.7 x 10~^*gcm~'^ (i.e., 
rtH = 3.4cm^'^ and tihc — 0.31 cm^'^). The other param- 
eters are set to the values inferred above: nno/nn = 1, 

"Ho" /"-He = 1, "-Hc+/"Ho = 0, kTns = 50 CV, i;^! = 

100 kms^^. The velocity is appropriate for i = 60° and 
d — 60 pc (Sect. ^^), which are the values used for the 
comparison with the observations (Fig. ||). 

In the panels with the fractional abundances in Fig. ^ 
one sees that these follow the pattern expected from the 
analytic solution of the structure of the nebula (Eq. |l8|) . 
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Fig. 5. Comparison of the observations with the predicted Balmer emission for an ionisation nebula. As Fig. but with 
the predicted Ha image shown in the top-right panel, and with modelled cross cuts overdrawn. The model reproduces 
the head of the nebula quite well, but underpredicts the brightness in the tail. This may be because in the tail the 
assumption of case A for the line emission is inappropriate (Sect. ^.3.4 ). 



Only behind the neutron star there are deviations, due to 
the expansion of the nebula. The temperature also follows 
the predictions ahead of the neutron star, but behind it 
the compression wave caused by the expansion is visible. 
Due to the expansion, the density is reduced by roughly 
a factor four far behind the neutron star. The velocities 
associated with the expansion are predominantly perpen- 
dicular to the axis of symmetry, as expected. The maxi- 
mum radially outward velocity is 6.5 kms~^ at 2 x 10^^ cm 
behind the neutron star. At large distances, this reduces 
to 3 kms^^ (the reduction is due to cooling; without cool- 
ing, the maximum outward velocity is about lOkms^^ 
and remains this high due to continued photo-ionisation) . 
No shocks form, since the velocity remains well below the 
sound speed everywhere: (u/cs)max = 0.3. 

The predicted Ha, H/3 and H7 fluxes are compared 
with the observations in Fig. |. As one can see, the agree- 
ment near the head of the nebula is very good, but near 
the tail the brightness is underpredicted. 

5.3.3. Fiddling with the parameters 

To see how sensitive our results are to the choice of pa- 
rameters, we also ran models for different choices, each 
time allowing ourselves the freedom to adjust pism and 
/bb such that the match to the observations was as good 
as possible. We treat the different parameters in turn. 

First, we consider the properties of the interstellar 
medium. If the undisturbed temperature were lower, all 
temperatures would be somewhat lower and hence the 
emission less efficient. However, a small change in /bb will 



compensate for this: for a medium temperature of 100 K, 
/bb = 2 will lead to a model that is hard to distinguish 
from our base model. The medium may also not be com- 
pletely neutral, i.e., ^0 = n-no/nn < 1. This again does 
not lead to major differences, as long as one chooses a 
density larger by a factor ~ i-e., keeps njjo approx- 
imately constant. Even for ^0 = 0.5, the results are not 
much different, except that the tail of the nebula becomes 
slightly narrower because of the "dead mass" that has to 
be dragged along in the expansion. 

Second, we varied the temperature of the neutron star. 
This will lead to a change in the number of ionising pho- 
tons as well as a change in their mean energy e. We 
find, however, that we can compensate with a suitable 
change in f^h- one requires /bb — 2.1 and 1.4 for tem- 
peratures fcTns — 40 and 60 eV, respectively. The results 
look very similar, except that for lower (higher) temper- 
ature, the tail is somewhat longer (shorter) and narrower 
(wider), and the head is very slightly dimmer (brighter). 
Compensation for the change in brightness requires a 
change in pi^m of ~5%. 

Third, we considered different inclination. We found 
that even for inclinations as low as 30°, the main effect 
is on the width of the tail, which becomes narrower; the 
length of the tail is hardly affected, while the brightness 
of the head of the nebula changes only slightly (brighter 
for lower inclination) . The somewhat counter-intuitive re- 
sult for the width of the tail is a consequence of the fact 
that the expansion velocity is inversely proportional to 
Vns (Eq. 29), and, therewith proportional to sini (since 

Vns = Vns.±/ Siui = d^ns/smi). 
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Finally, we considered the distance. For larger dis- 
tance, less dense interstellar medium suffices to reproduce 
the brightness of the head, as expected from Eq. ^ (e.g., 
Pigjjj = 5.4 X lO^^^'gcm^^ for d — 100 pc). One also finds 
that the tail of the nebula becomes narrower with increas- 
ing distance, as a consequence of the higher space velocity 
of the neutron star implied by the measured proper mo- 
tion. Indeed, the width decreases roughly inversely pro- 
portional with distance; for d — 100 pc, it is substantially 
narrower than observed. 

From our attempts, we conclude that the model re- 
sults are remarkably insensitive to fiddling with the pa- 
rameters. This is good in the sense that apparently the 
observed shape of the head of the nebula follows quite nat- 
urally from the model, but bad in the sense that, when one 
matches the brightness of the head, just as naturally one 
predicts a tail to the nebula that is too faint. We discuss 
below whether this might be due to one of the assumptions 
underlying our model being wrong. 

5.3.4. Verifying the assumptions 

For models that reproduce the brightness of the head of 
the nebula, the brightness in the tail is lower than ob- 
served. Thus, in its present form, the model is incorrect. 
This may relate to one of the assumptions, which per- 
haps lead us to overestimate the flux near the head of the 
nebula, or underestimate the flux in the tail. It may be 
that we miss an important thermal or emission process, 
but we argued above that photo-ionisation of Hydrogen 
and Helium, and collisional excitation and ionisation of 
Hydrogen dominate all other processes. 

The one assumption we have not yet verified, is that 
of case A for the line emission, i.e., whether the medium 
is sufficiently optically thin to Lyman lines that, e.g., any 
Ly/3 line produced by collisional excitation escapes the 
nebula before being converted in an Ha photon (plus two 
photons in the 2s-ls continuum). With temperatures and 
densities from our model in hand, we can verify whether 
case A actually holds. 

For the Lyman lines, the line-centre cross-section is 
given in Eq. |l^. At the position of maximum brightness 
in the nebula, 4.5 x 10^^ cm (0'.'5) in front of the neutron 
star, our model has nn" = 1.4 cm^"^ and T = 5.2 x 10^ K. 
Thus, over the size of the head of the nebula (1" or 
~ 10"'^^ cm), the line-centre optical depth for Ly(3 photons 
is Ti.y0,o = 5.8. Averaged over the Doppler profile, the 
probability that the photon escapes can be estimated by 
Peso — 1-72/(to -I- 1.72) = 0.23 (here, a spherical geome- 
try is assumed; Osterbrock 1989| , p. 101). Since every time 
a Ly/3 photon is absorbed, there is a probability of only 
Pconv.i = 0.118 that it is converted to Ha (instead of being 
scattered) , the total probability that an Ha photon is pro- 
duced is Pconv.tot ^ Pconv,l/(l " (1 " Pconv,l)(l ^ Pose)) = 

0.37. This is still substantially smaller than unity and 
shows that case A at least is a better choice than case B. 
Some of the photons may still be absorbed close to the 



head of the nebula, but many may also escape altogether, 
since in the hot central regions the Doppler profile is 
substantially broader than in the cooler surroundings. 
Obviously, for Ly7 and Ly5, the optical depths are lower 
and hence the case A assumption is better. 

The situation is somewhat different, however, in the 
tail of the nebula, where the neutral Hydrogen density is 
higher and the temperature lower. For instance, at z = 
1.8 X 10^^ cm (20" behind the neutron star; slit s), maxi- 
mum emission is at p = 5.0 x 10^^ cm (5'.'6 away from the 
nebular axis), and at this position one has njjo = 2.7 cm~^ 
and T = 2.0 x 10"* K. Over the relevant length scale at this 
position, ~3 X 10^^ cm, the optical depth is Ti^ypfl — 54. 
Hence, posc = 0.03 and Pconv.tot = 0.8. 

From the above, we see that case A is not quite appro- 
priate. If one assumes case A, therefore, the emissivities 
are underestimated, and therewith the required density 
overestimated. For the head of the nebula, where the gas 
is hot and highly ionised, the effect will not be very large, 
but in the tail, where the gas is relatively cool and less 
ionised, case A is clearly inappropriate. To treat this cor- 
rectly, one needs to do full radiation transport, which is 
beyond the scope of the present paper. To get an idea of 
the effect, however, we can make the opposite assumption 
for the tail, namely that the medium is optically thick to 
Lyman lines everywhere (i.e., case B). For this case, the 
emission in the tail would be increased by a factor ~ 2, 
which brings it into much closer agreement with the ob- 
servations. We intend to pursue this in more detail in a 
future publication. 

One other effect we have ignored so far, is con tin- 
uum pumping by the ionising source ( iFerland 19991) , m 
which neutrals absorb continuum photons at Lyman tran- 
sitions and emit Balmer, Paschen, and other higher-level 
lines. For Ha, we estimate a maximum photon rate of 
10~^ s~^ cm~^ from this process; for this estimate, we cal- 
culated the neutron star flux at Ly/? from the observed 
optical flux, and assumed that all photons within a ther- 
mal line width of 0.1 A (appropriate for T = 5 x 10^ K) 
were absorbed and converted to Ha (i.e., case B). Even 
this maximum rate is well below that produced by col- 
lisional excitation near the head of the nebula. For the 
Helium lines, however, continuum pumping will be impor- 
tant, as the nebular emission is due to recombination and 
therefore very weak. Given the rather large size of the 
region in which Helium is ionised, however, the emission 
will be rather diffuse, and hence it is no surprise we do 
not detect it. 



5.3.5. Model inferences 

We found that there were two main free parameters in 
our model, the interstellar medium density pism and the 
factor /bb, which measures the deviation of the spectrum 
from a black-body. For the density, the value depends on 
the fraction = nno/nn; for a distance of 60 pc, we 
found that pism — 7.7 x 10^^^£,q^ gcm^"^. Likely, the real 
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value is somewhat lower, since in assuming case A we ig- 
nored conversion of Ly/3 to Ha photons, which we found 
was not quite correct. Full radiative transfer calculations 
are required to infer a more precise value. Pending those, 
we round down the required Hydrogen number density to 
riH — 3^1^^ cm~^. 

The factor /bb depends mainly on the neutron-star 
temperature. This is because these two together set the 
rate and energy of the ionising photons, the combina- 
tion of which sets the total energy input, which is well 
constrained by the size of the nebula. Roughly, one has 
Ein (X nxaion(e ~ x) oc /bb^*^'^^. To keep ii^in constant, 
thus requires /bb oc T""'^^. Our model results confirm 
this expectation; we find that roughly /bb = I.TT^q^. 

6. Ramifications 

RX J1856.5— 3754 is the nearest neutron star and the 



brigh test soft X-ray source (in this category; Treves et al 
2000|). As such, it is commanding great interest, in particu- 
lar from spectroscopic X-ray missions. The available data 
inform us that RX J1856.5— 3754 is nearby (60 pc) and 
warm (50 eV), and that possibly it was formed in a nearby 
stellar association about a million years ago. Beyond this, 
however, we are very much in the dark about its nature. 

The source is puzzling for two reasons in particular. 
First, using the temperature inferred from the X-ray, ul- 
traviolet and optical energy distribution, the distance de- 
rived from the measured parallax and observed flux imply 
a radius of only ~ 7 km. This is too small for any equa- 
tion of state dPons et al. 200 1| ). Second, RX J1856.5-3754 
does not show any pulsations. The lack of pulsations in the 
radio could be due to the beam not intersecting our line- 



of-sight. The lack of X-ray pulsations (Pons et al. 2001 



Burwitz et al. 2001), however, is puzzling, since all well- 
studied cooling neutron stars do show significant pulsa- 
tions. 

In this paper, we report the discovery of a cometary 
nebula around RX J1856.5— 3754, shining in Ha and 
H/3. We have considered two models to account for the 
cometary nebula: a pulsar bow shock model and a photo- 
ionisation model. Both models lead to new, though differ- 
ent, diagnostics that bear upon the origin of this interest- 
ing object. 

We first summarise the diagnostics below (Sect. |6.l|) 
and then describe what future observatio nal t ests could 
distinguish between the two models (Sect. |6^ ). We con- 
clude with a brief discussion about the puzzling absence 
of pulsations (Sect. |6.3| ). 

6.1. Constraints Provided from Observations of the 
Nebula 

In both models, we can show that the pulsar is moving 
through a low density medium, nno — 0.8 cm~^ (bow- 
shock model) or n^o ~ 3cm~'^ (photo-ionisation model). 
Such densities are typical of the warm neutral medium 
(WNM) and perhaps even the warm ionised medium 



(WIM). These two phases of the interstellar medium are 
expected to pervade about half of interstellar space. Thus, 
it is not unusual to find RX J1856.5— 3754 embedded in 
one of these. 

The two models do not directly constrain the neutral 
fraction of the ambient gas, ^o- The WNM is expected to 
be essentially neutral, while for the WIM one finds <^o ^ 
0.5 (at least in those few cases where sufficient data exist 
to allow estimates; see Sect. 

Below we discuss additional diagnostics that are pecu- 
liar to each of these two models. 

6.1.1. Pulsar bow shock 

"Balmer-dominated" nebulae have been seen around a few 
pulsars. The pulsar wind combined with the supersonic 
motion of the pulsar in an ambient (partially) neutral 
medium results in a bow shock. The Balmer lines arise 
from the coUisional excitation of ambient neutral atoms 
as they drift into the post-shock region. 

In the framework of the pulsar bow-shock model, 
RX J1856.5— 3754 is a young off-beam pulsar with an en- 
ergetic pulsar wind powered by rotational energy loss E ~ 
6 X lO'^^^fj'^dgQ ergs^^. Combining this E with the inferred 
age of ~ 10^ yr, we estimate the spin period (P) and the 
strength of the magnetic dipole field {B): P ~ 3x^/^^y^ s 

and B ~ W^xil'^ G; here x = 1 - {Po/P? with Pq be- 
ing the spin period at birth. Since x 5; 1 ^-J^d S^o < 1, the 
estimates are upper limits; they are consistent with values 
for P and B expected from an ordinary million-year old 
pulsar. 

How secure are these conclusions? At the speed rele- 
vant to this source - lOOkms"^ - there currently exist no 
rigorous calculations of Balmer-dominated shocks. In ad- 
dition, we have not taken into account the pre-ionisation 
of the ambient gas ahead of the neutron star by the radi- 
ation from the neutron star and from the shocked gas, we 
have not properly modelled the bow shock geometry, and 
we have ignored issues of radiative transport. Thus, the 
constraints we derive must be viewed with some caution. 
With some effort, however, the full modelling effort can 
be undertaken (J. Raymond, pers. comm.). 

6.1.2. Photo ionisation nebula 

We also investigated the possibility that the nebular emis- 
sion is caused by photo-ionisation and heating of the ambi- 
ent gas by the extreme- ultraviolet radiation of the neutron 
star. Such nebulae were predicted to exist around rapidly 
moving, hot neutron stars by Blaes et al. (1995). They 



should have cometary shape and be rather different from 
normal ionisation nebulae, in that both recombination and 
metal-line cooling are irrelevant. 

We made analytical estimates which showed that the 
nebula around RX J1856.5— 3754 could indeed be due to 
photo-ionisation. Since the input atomic physics for the 
photo-ionisation model is well known and relatively sim- 
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pie, we proceeded by building a detailed model. We find 
that the model can reproduce the observations in some de- 
tail. There are discrepancies as well, but we argued that 
those could be related to the fact that for the line emission 
neither case A nor case B is applicable, and hence one has 
to do full radiative transport. 

We concluded that an ionisation nebula could repro- 
duce the observations, provided the ambient density is 
about S^Q^^cm^^ and that the extreme- ultraviolet flux 
(which is responsible for most of the photo-ionisation) is 
somewhat higher, by a factor ~1.7, than expected based 
on the black-body spectrum that best fits the optical- 
ultraviolet and X-ray fluxes. Since some deviations of the 
spectrum from the blackbody distribution are expected, 
and since deviations are indeed seen at X-ray wavelengths, 
the required excess at extreme ultra-violet wavelengths 
does not seem unreasonable. Indeed, in the context of this 
model, the nebula provides a measurement of the extreme 
ultraviolet flux, something which is not possible otherwise 
because of interstellar extinction. 

Proceeding further, in the photo-ionisation framework 
(but not in the bow shock model), it is plausible that 
the ambient gas will accrete onto the neutron star (if 
not stopped by a pulsar wind or flung off by a rotat- 
ing magnetic field). Assuming the Bondi-Hoylc mecha- 
nism, we find an accretion radius of i?acc = 2GMns/u^j, = 
3.7 X W^d cm. This leads to an accretion rate of 



accPismi'ns 
-0.5 , 



3 X lO^dgQ' gs (where we used 



nu = 3dgQ'^ '^cm and X — 0.73). The resulting accre- 

6 X lO^^detf-^ergs-^ In 



tion power is ~ GM^gM / R^, 
contrast, the bolometric luminosity of RX J1856.5— 3754 
is ^4 X 10^^d|o ergs"^ 

Thus accretion, even if takes place, is energetically in- 
significant. RX J1856.5— 3754 is hot not because of accre- 
tion. However, accretion will also introduce cosmic com- 
position material to the surface of the neutron star. If the 
metals settle out rapidly, as is expected for accretion at 
levels too low to contribute significantly to the luminosity 



( [Brown et al. 1998| ; [Bildsten et al. 1992| ), this would lead 
to an atmosphere of almost pure Hydrogen. This would 



be inconsistent with the observations (Pons et al. 2001 
Burwitz et al. 20011). 



6.1.3. A combination? 

If the nebula is indeed formed by photo-ionisation then 
any E associated with rotational energy loss should be 
small enough that the standoff radius, ra, is less than the 
observed size of the nebula. Adopting the density required 
by the ionisation model, rin — 3cm~^, this leads to an 
upper limit, E 2 x 10'^^dg5^ ergs~^. As before, combin- 
ing knowledge of E with the age of 10^ yr inferred from 
an origin in Upper Scorpius, we derive P ^ ^x^^"^ s and 
B ^ 7 X lO^^xG. Unfortunately these latter constraints 
are weak, since < x < 1. 

Conversely, if the factor /bb is significantly larger than 
unity, /bb <^ 4, then the bow shock model is ruled out. 



6.2. Observational Tests 

The critical difference between the pulsar bow shock 
model and the photo-ionisation model is in the velocity 
field. In the former, we expect to see more than half 



(Raymond 1991) of the Ha emission to have velocity 
widths comparable to the shock speed of lOOkms"^ both 
because of the high post-shock temperature and because 
of large bulk motion. In contrast, in the photo-ionisation 
model the thermal velocities are moderate, ^40kms^^, 
and the bulk motion even smaller, < 10 km s~^ (see Fig. ||). 
Thus the simplest test is to obtain higher spectral resolu- 
tion data. If only narrow lines are seen, then the photo- 
ionisation model is vindicated. With adequate spatial res- 
olution, one can map the velocity field, for which again 
the expectations are clearly different in the two models. 

A high spatial resolution image of the nebula would 
be extremely useful. First, independent of which model 
is correct, such an observation can be used to derive the 
direction of the velocity vector of RX J1856.5— 3754 (be- 
cause of the intrinsic cylindrical symmetry of both mod- 
els). Next, it can be used to see whether the nebula is thin 
with a sharp edge, as expected in the bow shock model, or 
whether its emission gradually peters out, as expected in 
the photo-ionisation model. A pulsar wind would evacuate 
gas up to its stand-off radius and thus the most interest- 
ing possible result of a high resolution image would be the 
discovery of a cavity around RX J1856.5— 3754, contained 
within the nebula. If so, the photo-ionisation model would 
be vindicated and we would still obtain a direct measure 
of 



6.3. Lack of X-ray Pulsations 

From the proper motion, it seems likely that 
RX J1856.5-3754 originated in the Upper Sco OB 
association and that it is a a million-year old cooling neu- 
tron star. Indeed, for the bow-shock model, the inferred 
value of an E is typical of a million year old pulsar. If 
so, why does RX J1856.5— 3754 not exhibit pulsations? 
A radio beam might not intersect our line-of-sight, but 
all well-studied cooling neutron stars show significant 
pulsations at X-ray wavelengths. 

The pulsation in cooling radiation arises because mag- 
netic fields influence the conductivity in the crust and 
thereby lead to temperature variation at the surface. A 
hotter thermal component (in addition to any magneto- 
spheric emission) is seen in energetic pulsars and usually 
attributed to a polar cap heated by pulsar activity. 

Among the middle-aged pulsars, we note with inter- 
est that the one with the weakest field, PSR B1055— 52 
(P = 197 ms, B - 1.1 X 10^2 Q and ~ 3 X lO^^erg 
s^^), is also the one whose flux show s weakest modula- 



tion, ^ 10% (lOgelman fc Finley 1993]) . Could it be that 
RX J1856. 5-3754 is a slightly older analogue of PSR 
B1055— 52? The pulsed fraction in the cooling emission 
for RX J1856.5— 3754 could be reduced further by appeal- 
ing to an even weaker magnetic field strength, say to a few 
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10^^ G (in the bow-shock picture, this would imply that 
X ^ 1, i.e., that the pulsar was born with a period close 
to the current one). If so, it might have been born with 
properties like those of the 39.5-ms pulsar PSR B1951+32, 
which has _B ~ 5 x 10^^ G and is associated with the su- 



pernova remnant CTB 80 (Kulkarni et al. 1988). 

If the suggestion that RX J1856.5— 3754 has a mag- 
netic field of a few 10^^ G is correct, we expect the follow- 
ing: (1) more sensitive observations will reveal pulsations 
at the level of a few percent; (2) E, when measured, will 
be below 2 x lO^^ergs"^; (3) at high sensitivity, a weak 
hot component, ^ 1% of the total emitted flux, may be 
identified. All these expectations can be tested with deep 
Chandra and XMM observations. 

Most pulsars that we observe, however, appear to 
be born with magnetic field strengths upward of 10^^ G. 
In this respect, the a priori probability seems low that 
RX J1856.5— 3754 has such a weak field. It may be that, 
after all, RX J1856.5— 3754 did not originate in Upper Sco, 
and is old, lost its magnetic field due to accretion, and 
perhaps was recently reheated by passing through a dense 
molecular cloud. Or perhaps RX J1856.5— 3754 was just 
born with a very low magnetic field strength, or lost its 
field in a binary with a complicated history. These possi- 
bilities seem even less plausible. Thus, if our chain of logic, 
which admittedly is long and has many weak links, is cor- 
rect, RX J1856.5— 3754 is a middle-aged relatively weakly 
magnetised pulsar whose beam does not intersect our line 
of sight. 

It is ironic that the two brightest nearby neutron stars, 
RX J1856.5-3754 and RX J0720.4-3125, may well rep- 
resent the extreme ends of the neutron star magnetic fleld 
distribution, one a weak field neutron star and another a 
magnetar. 
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